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Abstract 

We consider the slave-fermion Schwinger-boson decomposition of an effective 
model obtained through a systematic low-energy reduction of the three-band 
Hubbard Hamiltonian. The model includes a three-site term t" similar to 
that obtained in the large-t7 limit of the Hubbard model but of opposite sign 
for realistic or large 0-0 hopping. For parameters close to the most realistic 
ones for the cuprates, the mean-field solution exhibits d + s superconductivity 
(predominantly d x i_ y i) with a dependence on doping x very similar to the 
experimentally observed. We also obtained incommensurate peaks at wave 
vectors ~ 7r(l, 1 ± 2a?) in the spin structure factor, which also agree with 
experiment. 

PACS 74.20: Theory of superconductivity. 

75. 10. J: Heisenberg and other quantized spin models. 
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After the original derivation of the t — J model by Zhang and Rice there has been 
much discussion about its validity 0-§|. At present it seems clear that while the low-energy 
eigenstates contain ~ 98% of Zhang-Rice local singlets ||] supporting the reduction of the 
Hilbert space, other terms in addition to the nearest-neighbor hopping t\ and exchange J 
should be included to accurately describe the low-energy physics of the original three-band 
Hubbard model |U|[7]]. In particular, numerical comparison of the energy levels of a CU4O8 
cluster has shown that it is necessary to include the three-site correlated-hopping term t" 
(see Eq. (1)) to obtain a reasonable fit of the low-energy spectrum [|]]. 

While in the one-band Hubbard model t" < (with the notation of Eq. (1)), in the 
three-band Hubbard model t" changes sign with increasing direct 0-0 hopping t pp and 
becomes positive for realistic or large t pp |5]||. Since most of the terms in t" have the 
effect of rotating 90° a nearest-neighbor singlet, one might expect that a negative (positive) 
t" favors s (<i)-wave pairing in the ground state. This is in agreement with Monte Carlo 
calculations which show indications of s-wave (d-wave ||10|| ) superconductivity in the 
three-band Hubbard model for t pp = (sizeable tpp). In addition, exact diagonalization of 
a cluster of 4 x 4 unit cells has shown important indications of a sudden large increase in 
the G?-wave superconducting correlations and off-diagonal long-range order for t" > t", where 
t" c jt\ ~ 0.1 ||. Simultaneously the wave function becomes similar to a simple short-range 
(d + s)-wave superconducting resonance-valence-bond (RVB) wave function. The physics 
behind this is simple: the RVB wave function takes advantage of both, the short-range 
magnetic correlations induced by J and the hopping of nearest-neighbor singlets t". 

However, the above mentioned studies suffer from finite-size effects and either the tem- 
perature or the doping are not realistic. In particular, in the 4x4 cluster, adding only two 
pairs of holes the system is already in the overdoped region of real materials x = 0.25, and 
a value of t" nearly three times larger than the one obtained from the mapping procedure 
was necessary to enter the RVB superconducting region [[|. 

In this work we use a slave-fermion Schwinger-boson treatment |TI]-|rjJ of the effective 
t — t" — J model. While being approximate, this method does not have the above mentioned 
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shortcomings. In the undoped system the effective model reduces to the Heisenberg model, 
whose ground-state properties are well accounted for by the Schwinger-boson approximation 



i~T| , |T5l| . For the doped system, we calculate the doping-dependence ot both the supercon- 
ducting gap and the wave vector of the observed incommensurate peaks in the magnetic 
fluctuations |TjJ, obtaining a good overall agreement with experiment. While previous nu- 



merical work in the one-band Hubbard JT^| and generalized t — J model fl8| obtained a 



deviation with doping x of the maximum position in the spin structure factor from (7r, ti) 
towards (n, 0), to our knowledge no theory so far has been able to explain the fact that the 
deviation amount is very near to 2tcx |16| . 



The model can be written as 



H=-J2t ir clc i<7 + t" £ c5, a c i ^-2S r S j ) + ^(S r S i -^ i n j ). (1) 



ira i,j¥ : j'> cr ~~ V a 

The first term contains hopping to first, second and third nearest neighbors with parameters 
ti, t 2 , t 3 respectively. In the second, correlated hopping term, j and f label the four 
nearest neighbors of site i. The last term is the exchange term. The operator c\ a creates an 
electron on a vacuum where all sites carry a Zhang-Rice singlet HU and double occupancy 
is forbidden (of^n, = 0). The parameters ti, t" and J are taken within the (narrow) band of 
possible values derived from a realistic three-band Hubbard model ]5|||||. We take t\ — 1 
as the unit of energy and J = 0.4, t" = 0.05. Instead, t 2 and t 3 are included to account 
for effective second and third nearest-neighbor hoppings which are dynamically generated 
in the bare t — J model, and which are not taken into account in our mean- field treatment. 
With t 2 = 0.25, t 3 = 0.18, the quasiparticle dispersion of the bare t — J model calculated 
using the self-consistent Born approximation [19j or Monte Carlo ]2(J is well reproduced by 



our approximation. 

The constraint of no double occupancy is treated using the slave-fermion Schwinger-boson 
decomposition, q ct = bi a fj , c\ a = fib\ a which, together with the site-occupation constraint 
Sct&LAct + fifi = 1, is a faithful representation of the original Fermi algebra. In this case 
the spins can be written in terms of the Schwinger bosons as Si = \b\ a a a f3bip, while the 
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slave fermion fa accounts for the charge degrees of freedom. We prefer this decomposition 
since the Schwinger bosons are known to produce a very good description of the long-range 
antiferromagnetic order at half filling and the incommensurate magnetic fluctuations upon 



doping [|13|-|15|. The proper account of these magnetic fluctuations are essential to describe 
the low-energy physics of the model, both at and near half filling. After the Hamiltonian 
is written in terms of the new operators, the partition function is expressed as a functional 
integral over coherent Bose and Fermi states. Furthermore, the different terms of Eq. (1) 
are decoupled by means of Hubbard-Stratonovich transformations according to the following 
scheme (indicated by a vertical bar): 

H t — -EwjJ i bib raj H t „ — t" J2 01 f] \frfi - fjfr) 1 4- l K 



ira i,j¥"j' 

Hj^jEt 4 I 4 = "4 I 4 - (1 " *?) " | E {flfj I fifi 1 4 1 4) (2) 

<ij> <ij> 

The rationals behind these decompositions are as follows: For the hopping terms we just 
decoupled bosons from fermions. In Hj we took into account the pairing effects over the 
holons / of breaking magnetic bonds, and we wrote the Schwinger-boson spin-spin inter- 
action in terms of the two SU(2) invariants Bjj = ^a-bl^bju, Aij = J2a ^b ia bj^. These 
operators describe the magnetic fluctuations in the ferromagnetic and antiferromagnetic 
channels respectively, and allow for a proper treatment of the incommensurate order upon 
doping. Finally, H t n is conveniently written as boson and fermion pair hopping, in terms 
of the antiferromagnetic singlet and holon pairing operator. This decoupling highlights 
the connection between superconductivity and short-range magnetic fluctuations observed 
in numerical studies of the Hamiltonian ||. 

According to the above decomposition of the interaction terms, in a saddle-point evalu- 
ation of the partition function, the decoupling Hubbard-Stratonovich fields become the fol- 
lowing self-consistent order parameters (77 connects nearest neighbors): 1) antiferromagnetic 
order: A v = (5Z<j °"&i<r&i+»7,er) ) 2) ferromagnetic order: B v = bl a bi +rit(T ) , 3) superconduct- 
ing order: D v = (fifi+ v ), an d 4) nearest-neighbor hopping: F v = (f-f i+r) ). In k-space the 
saddle-point Hamiltonian can be written as: 



k,cr 
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The subscript k refers to the usual Fourier transform. However, we have defined b[ 
"^Sk^ko- e lk<Tri with k CT = k + crQ/2. The extra phase crQ.r/2 is required to obtain 
the gapless (Goldstone) modes associated to the long-range magnetic order at the proper 
k = 0, ±Q points. In Eq. (||) we have defined: 

A k = A + E trjF^ + -B^ R k = -A k + 2t" £ F^A^' 
v 4 4 nfrf 

-2t"L> k A V D* V + i {At" + J){\D X \ 2 A X sin k x + \D y \ 2 A y sin k x ) , 



T k = E - 2i (*" + j)(\ A x\ 2 Dx sin ^ + l^l 2 ^ sin k y 

v 



(4) 



The Hamiltonian (^]), quadratic in boson and fermion operators, can be diagonalized in the 
standard way. The order parameters are determined selfconsistently using a mesh of 200x200 
points for the integrals in reciprocal space. Since our mean-field theory is not variational 
because it unphysically enlarges the configuration space, we cannot rely only on the energy- 
minimization criterion to seek for solutions of the consistency equations. Consequently, 
based on previous numerical studies of the t — J model |I7| , |18| we sought for solutions 
corresponding to magnetic wavevectors (n ± Q, tt) (which in some parameter region can in 
fact have slightly higher energy than the corresponding to a wavevectors along the diagonal). 
For these solutions, the values of the bosonic order parameters as a function of doping x are 
shown at the left of Fig. 1. For i ^ 0, the Neel order turns continuously into a spiral spin- 
density wave of increasing pitch, as a consequence of both the finite value of the ferromagnetic 
order parameter B x and the decrease of the antiferromagnetic parameters A x , A y (mainly in 
the direction of the spiral x). As expected, the dynamical susceptibility x(q, lj) is peaked 



at wave vectors (n ±Q,tt). The dependence of Q on doping is very approximately linear, 
as represented at the right of Fig. 1. The ratio r = Q/(2ir) ~ 1.2, and is weakly dependent 
on the parameters. It decreases for larger J. The experimentally observed incommensurate 
magnetic fluctuations (with r ~ 1) |TJ| are consistent with microdomains of spiral density 
waves of both possible directions (x and y). 

The superconducting order parameter is given by 

( c ikT c -k| - c k| c -kt) = 2 ( A * cos k x + A y cos k y ) , = A V D V . (5) 

The values of A^ are represented in Fig. 2. To show the effect of t" we also include in 
this figure the corresponding result for t" = and t" = —0.05. In agreement with the 
experimental dependence of T c with doping, the absolute value of the order parameter is 
maximum for a doping slightly larger than x = 0.15. As a consequence of the symmetry 
breaking induced by the spiral spin arrangement, |A X | ^ \A y \ and the order parameter has 
always mixed s + d character. For t" > 0, A x and A y have opposite signs, and the d- 
wave component dominates. In agreement with previous numerical studies 0, a positive t" 
increases considerably the magnitude of the order parameter. The present results indicate 
that a small, realistic t" leads to the proper order of magnitude of the superconducting 
gap (~ (J + At") A y ~ 0.005eV for optimum doping). Instead, a negative value of t" (in 
the notation of Eq. (1)) suppresses superconductivity. In particular, for t" < —0.04 the 
signs of A x and A y are the same and, therefore, the s-wave component dominates the small 
superconducting signals. Since important direct 0-0 hopping in the three-band Hubbard 
model leads to positive t", while t" = —J/4 in the large-f/ limit of the Hubbard model, 
the present results agree with Monte Carlo results which suggested s-wave (c?-wave [[[(J ) 
superconductivity in the three- band Hubbard model for t pp = (sizeable t pp ) and the 
absence of strong superconducting signals in the one-band Hubbard model. 

In conclusion, using a slave-fermion Schwinger-boson mean-field approximation applied 
to an effective model for the cuprates, we obtain a superconducting state of mixed (d + s)- 
wave character with approximately 75% of <i-wave component. The order of magnitude 
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of the gap and its dependence with doping agree qualitatively with experiment. More- 
over, he symmetry of the order parameter agrees with different experiments which show the 
existence of nodes, opposite signs along x and y directions, but also important tunneling 
currents along the c axis in junctions with conventional s-wave superconductors [^TJ. In fact, 
Annet et al. have shown that a real combination of s and c?-wave components is consistent 



with most experiments |22[. However, it is necessary to break the tetragonal symmetry to 
avoid a splitting of T c . We also obtain a spiral spin-density wave which provides this sym- 
metry breaking and is consistent with the incommensurate magnetic fluctuations observed 



in superconducting L^-xSr^CuC^ [[L6[]. The spiral spin density wave can be stabilized by 
the effective negative next-nearest-neighbor hopping which appears in the holon dispersion 



relation |13|Jl^] , in agreement with numerical studies JT^JTMI • Another possible explanation 
of the incommensurate fluctuations in terms of dynamical stripe phases has been suggested 
23| , |24|| . A static version of this phase consists of domain walls parallel to the x (or y) 



direction containing one hole every two Cu sites, and is consistent with neutron-diffraction 
experiments in Lai. 6 _ x Ndo.4Sr a; Cu04 with x ~ 1/8 p3[] . However, domain- wall solutions of 
Hubbard or t — J models are always insulating P5| , P^] , and as stated by Schulz |]25[| , motion 
of domain walls is unlikely because they are easily pinned by impurities and electron-lattice 
interactions. Although further study is needed to establish the origin of the observed in- 
commensurate magnetic fluctuations, we find it encouraging that our approach is consistent 
with them and with the essential features of the superconducting state for the same set of 
realistic parameters. 

Helpful discussions with J.J. Vicente are gratefully acknowledged. Three of us 
(CDB,HAC,AAA) want to thank the hospitality of the International Centre for Theoreti- 
cal Physics, where part of this work was done. CDB and LOM are supported by Consejo 
Nacional de Investigaciones Cientificas y Tecnicas (CONICET), Argentina. AAA and HAC 
are partially supported by CONICET. 



7 



REFERENCES 

[1] F.C. Zhang and T.M. Rice, Phys. Rev. B 37 (1988) 3759. 

[2] V.J. Emery and G. Reiter, Phys. Rev.B 38 (1988) 11938. 

[3] C.X. Chen, H.B. Schiittler and A.J. Fedro, Phys. Rev.B 41 (1990) 2581. 

[4] CD. Batista and A. A. Aligia, Phys. Rev. B 48 (1993) 4212; 49 (1994) 6436 (E). 

[5] A. A. Aligia, M.E. Simon, and CD. Batista, Phys. Rev. B 49 (1994) 13061. 

[6] J. Eroles, CD. Batista and A.A. Aligia, Physica C 261 (1996) 237. 

[7] L.F. Feiner, J.H. Jefferson, and R. Raimondi, Phys. Rev. B 53 (1996) 8751; references 
therein. 

[8] CD. Batista and A.A. Aligia, Physica C 264 (1996) 319; J. Low Temp. Phys. 105 
(1996) 591. 

[9] G. Dopf, A. Muramatsu, and W. Hanke, Phys. Rev. Lett. 68 (1992) 353. 
[10] K. Kuroki and H. Aoki, Phys. Rev. Lett. 76 (1996) 4401. 
[11] A. Auerbach and D. Arovas, Phys. Rev. Lett. 61 (1988) 617. 

[12] B. Shraiman and E.D. Siggia, Phys. Rev. Lett. 62 (1989) 1564; D.M. Frenkel and W. 
Hanke, Phys. Rev. B 42 (1990) 6711. 

[13] B. Chakravorty, N. Read, CL. Kane and P.A. Lee, Phys. Rev. B 42 (1990) 4819. 

[14] B. Normand and P.A. Lee, Phys. Rev. B 51 (1995) 15519. 

[15] H.A. Ceccatto, C.J Gazza, and A.E. Trumper, Phys. Rev. B 47 (1993) 12329. 

[16] S.W. Cheong, G. Aeppli, T.E. Mason, H. Mook, S. Hayden, P.C Canfield, Z. Fisk, K.N. 
Clausen, and J.L. Martinez, Phys. Rev. Lett. 67 (1991) 1791. 

[17] D. Duffy and A. Moreo, Phys. Rev. B 52 (1995) 15607. 



[18] R.J. Gooding, K.J.E. Vos, and P.W. Leung, Phys. Rev. B 49 (1994) 4119; ibid 50 (1994) 
12866. 

[19] Z. Liu and E. Manousakis, Phys, Rev. B 44 (1991) 2414. 

[20] E. Dagotto, A. Nazarenko, and M. Boninsegni, Phys. Rev. Lett. 73 (1994) 728. 

[21] A.G. Sun, D.A. Gajewski, M.B. Maple, and R.C. Dynes, Phys. Rev. Lett. 72 (1994) 
2267. 

[22] J. Annet, N. Goldenfeld, and A.J. Legget, to appear in Physical Properties of High 
Temperature Superconductors, Vol. 5, D.M. Ginsberg (ed.) (World Scientific, Singapore, 
1996). 

[23] J.M. Tranquada, B.J. Sternlieb, J.D. Axe, Y. Nakamura, and S. Uchida, Nature 375 
(1995) 561. 

[24] J. Zaanen, M.L. Horbach, and W. van Saarloos, Phys. Rev. B 53 (1996) 8671. 

[25] H.J. Schulz, Phys. Rev. Lett. 64 (1990) 1445. 

[26] D. Poilblanc and T.M. Rice, Phys. Rev. B 39 (1989) 9749. 

FIGURE CAPTIONS 

Fig.l: Left: Order parameters A V ,B X as a function of doping x for one of the solutions 
(that with By = 0) which minimize the free energy (the other solution is equivalent by a 
rotation of 90°). Right: position of the maxima of the dynamical susceptibility expressed as 
(n ± Q, 7r) as a function of doping. 

Fig. 2: Superconducting order parameter (see Eq. (|5|)) as a function of doping for three 
values of t" (from top to bottom, t" = 0.05, and -0.05). 
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